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Vortices in 5*0(3) x Z(2) simulations 

Andrei Alexandru and Richard W. Haymaker* a 

a Dept. of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana, 70803 USA 

We explore simulations on periodic lattices in the Tomboulis SO(3) x Z(2) formulation. The dynamical vari- 
ables are constrained. We propose an update algorithm that satisfies the constraints and is straightforward to 
implement. We show how boundary conditions put constraints on the configuration space. 



1. Introduction 

The center vortices seem to play an important 
role in disordering the Wislon loop. The mech- 
anism is explained in a number of papers us- 
ing an intuitive idea about vortices. We decided 
to investigate the vortices using Tomboulis for- 
malism [1,2]. The Tomboulis formalism is based 
on spliting the SU(N) degrees of freedom into 
SU{N)/Z(N) variables living on links and Z{N) 
variables living on plaquettes. The new variables 
are constrained and for the free boundary case 
the constraint amounts to a coincidence between 
the thick and thin monopoles. 

In trying to simulate the new variables we have 
to solve two problems. The first problem regards 
the boundary conditions: the Tomboulis analysis 
was done for a free boundary lattice. We need 
periodic boundary condition for our simulation 
and we have to understand how they influence the 
constraint on the variables. The second problem 
is to find a suitable definition for the constrain- 
ing set. The definition that Tomboulis derived 
was rather awkward to implement in a numerical 
investigation. 

2. Derivation 

From now on we are going to deal only with 
SU(2). We start with the Wilson partition func- 
tion for SU{2): 



[dU]e 



■E v t < u p) 
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Following Tomboulis we re- write it: 

Z= [[dU}'J2C[<rr)}e^ lTr(Up)la{P \ 
where 

C H = EII^)( Q W)IIW)- 

t£A p b 

The set: 

A={a:P^Z(2)} 

is the set of all Z(2) configurations defined on 
plaquettes and r](p) — sgnTr(U p ) is the sign of the 
plaquette. We see that the function C is defined 
as a summation over all configuration r that obey 
the constraint n6^( T (^)) = 1 of the function 
UpXr(p)(a{p)) where: 

Xi(±l) = 1, 
X-i(±l) = ±1, 

are the characters of the two irreducible represen- 
tations of Z(2). 

The partition function is then given by an in- 
tegration over the SO(3) degrees of freedom that 
live on the links and a summation over the Z{2) 
degrees of freedom that live on plaquettes. The 
integration over the 5*0(3) is unconstrained. The 
summation over the Z(2) variables seems to be 
unconstrained too but the C function will act as 
a constraint. For example for free boundary con- 
ditions the C function has the form: 

C[a] =Y[S(a(dc)). 

c 

To make these things clear we will introduce 
some notation. 



3. Notation 

For any two configurations a, (3 G A we define: 



(a,/3) =Y[xa(p)(0(p))- 



This bracket is giving some degree of superposi- 
tion between the two configurations. We see that 
Xa(b) = — 1 only when both a and b are — 1. Thus 
the bracket is going to be 1 if we have an even 
number of plaquettes that have both a(p) = — 1 
and (3(p) — —1 and is going to be —1 if we have 
an odd number of such plaquettes. 

Here we list some properties of the bracket: 



(a, (3) -- 


-- (/?, a), 


a, /?7) = 


= (a,/3)(a,7), 


(a,l) = 


= 1. 



where 1 is the configuration with 1 (p) = 1 for all 
plaquettes p. 

In the definition of C we have a summation 
over all configurations that have the property 
n b <5(r(<%)) = 1. Define: 

C = {TGA\l[6(T(db)) = l}, 



the set of all such configurations. One important 
thing to note about the set C is that it forms a 
group under the multiplication law: 

(a/?)(p) = a(p)/J(p). 

Now for any subgroup /C included in A (C is 
such a subgroup) define: 

K = {a G A\ < a, [3 >= 1 V/3 G £}. 

In other words the set ^ is made up from the ele- 
ments that have the bracket with all the elements 
in the original set K, equal with 1. This is some 
sort of dual set since K, = K, and \fC\ x \)C\ = \A\ 
(we denote with |/C| the number of elements in 
the set K). 

Using this notation we will show that the C is 
just some sort of delta function and thus acting 
as a constraint. 



4. The C function 

Using our notations we write: 

C[a]=£>,r>. 

tGC 

Now we remember that C is a group and using 
group summation invariance we can write: 



C[a] 



rec 



<tt ) = (a,To)'y]{a,T) = (a,T )C[a\. 



rec 



This is true for any to € C. Thus if we have at 
least one to G C with (a, To) = —1 then C[a] = 
—C[a] and thus C[a] = 0. 

Now if we don't have any element in C that has 
the bracket with a equal to — 1 that means that 
(a, t) = 1 for all t G C and thus a is a element of 
the dual set C. Moreover: 

C[a]=5>,r)=£l = |C|. 

rec rec 

Summing up, we have: 



C[a] 



\C\ a € C, 

a^C, 



and we see that the function C is nothing more 
than the characteristic function for the set C (up 
to a constant). Thus the role of the C function in 
the partition function is to constrain the summa- 
tion over the Z{2) degrees of freedom to the set 
C. 

5. C set properties 

In the free boundary conditions 
case Tomboulis[l] showed that any configuration 
that obeys the cubic constraint n c ^( <J (^ c )) = 1 
is allowed. That is a configuration a is allowed if 
and only if for any cube c in the lattice we have 
an even number of plaquettes p in the boundary 
of the cube that have a(p) = — 1 (i.e. a(dc) = 1). 
Now we investigate the allowed configuration set 
in the case of periodic boundary conditions. 

In the periodic boundary conditions case all 
configurations a G C obey the cubic constraint. 
To prove it assume that : 

a(dc) = -1 



for some particular cube c. Then we take r G C 
defined as: 



r{p) 



•1 



1 



p G dc, 
p g - dc. 



Then (r, a) = — 1 since there are an odd number 
of plaquettes that have both a = — 1 and r = — 1 
(this is due to the fact that all plaquettes on the 
faces of the cube have r = — 1 but only an odd 
number of those plaquettes have a = — 1 since 
aide) = —1). Thus we proved that if the config- 
uration a doesn't obey the cubic constraint it is 
not a member of the C set. 

We see that even in the periodic boundary con- 
ditions case the allowed configurations have to 
obey the cubic constrained. However, this is only 
a necessary condition and is not sufficient as for 
free boundary conditions. To prove this take the 
configuration: 



a(p) = 



-1 peS 12 , 

1 P <£ Sl2, 



that obeys the cubic constraint. Now choose: 



Tip) 



-1 PGPi 2 , 
1 V $ P12, 



where S\i is a co-plane and P\i is a plane in the 
12 direction. Now it is easy to see that r G C 
and (t, a) = — 1 since S12 and P12 have only one 
plaqucttc in common. This means that although 
the a configuration obeys the cubic constraint it 
is not a member of the C set. 

Summing up these two observations we see that 
contrary to the free boundary conditions case the 
set C (the set of allowed configurations) is only a 
subset of the set of all configurations that obey 
the cubic constraint. Thus the periodic bound- 
ary conditions impose further constraints on the 
allowed configurations set. 

6. Alternative definition for the set C 

The partition function (up to a constant) is: 

z=[[dU}J2^ lTr(Up)Hp) . 

arjEC 



We see that the set C has all the informations 
regarding the constraint. However, the defini- 
tion that we have for the set C is not suitable 
for numerical simulations. The set C is defined in 
terms of the set C which in turn is defined using 
a constraint. Even in the case of free boundary 
conditions where the set C is given by the cubic 
constraint the definition is not easy to implement 
numerically. This is the reason we need a different 
definition for the set C. 

To do this we employ the star transformation 
around a link b defined as: 



a m(p) 



-1 p G db, 
1 p g" db. 



These are Z(2) configurations that are members 
of the set C. Since the set C forms a group any 
product of such configurations is a member of the 
set. Now we define the set of all star transforma- 
tions: 

V = {a G A\a = Y[aQ b .}, 

i 

which is obviously a subset of the set C. Now, 
the set T> forms a group and we can generate 
it by starting with the identity (all plaquettes 
equal with 1) and then doing star transformations 
around various links b (i.e. flipping the signs of all 
six plaquettes around a link at once). If we em- 
ploy this transformations for our update we are 
guaranteed to stay in the set C and thus obeying 
the constraint. Our only problem is to determine 
whether or not we are sweeping through the whole 
set C. We were able to prove [3] that V — C and 
thus we have proved that by doing the star trans- 
formations we are covering the entire set C. 
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